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Abstract 

We estimate the contribution of inelastic nucleon excitations to the (e,e') inclusive 
cross section in the CEBAF kinematic range. Calculations are based upon parameter- 
izations of the nucleon structure functions measured at SLAC. Nuclear binding effects 
are included in a vector-scalar field theory, and are assumed have a minimal effect on 
the nucleon excitation spectrum. We find that for q ^ 1 GeV the elastic and inelastic 
nucleon contributions to the nuclear response functions are comparable, and can be 
separated, but with roughly a factor of two uncertainty in the latter from the extrapo- 
lation from data. In contrast, for q ^ 2 GeV this uncertainty is greatly reduced but the 
elastic nucleon contribution is heavily dominated by the inelastic nucleon background. 
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1 Introduction 



Inelastic scattering of electrons from nuclear targets has long been a tool for the study of 
nuclear structure. For example, the Coulomb response has been experimentally determined 
with a view to testing the Coulumb sum rule against models of nuclear structure. Most 
data and analyses have been obtained for three-momentum transfers q < 550 MeV/c,[0-||^ 
although some recent work has been done for g ~ 1 GeV/c.p] At the lower q and energy 
transfers u {< q) the nuclear response involves almost entirely nucleons without internal 
excitation: the A-resonance is excited at higher u, but contributes mainly to the transverse 
response, and therefore has less effect on the Coulomb response. 

Interest in extending the study of the nuclear (e, e') response to higher q has developed 
in recent years, especially with the advent of the Continuous Electron Beam Accelerator 
Facility (CEBAF), with beam energies of ~ 4 GeV or more. This has led to the study of 
relativistic effects in the nuclear response which become important for g ^ 1 GeV/c.[3-||l6| 
But for such momentum transfers the probability of exciting internal states of the nucleons 
becomes increasingly important. It is the purpose of this paper to make an estimate of the 
magnitude of the contribution of internal excited states of the nucleon to the nuclear response 
based on measured values of the inelastic nucleon structure functions and a Fermi gas model 
of the nuclear target. The contribution of inelastically excited nucleons may be considered 
the background to the nuclear response function with non-excited (elastic) nucleons; it is 
the latter which is usually compared to models of nuclear structure. Our estimates are for 
1 < g < 4 GeV/c, with the appropriate ranges of u. 

Electroexcitation of nucleons has been studied at the Stanford Linear Accelerator (SLAG) 
and elsewhere in this momentum range, but more recently at much higher momenta. We 
make use of systematic fits|T^, |T^ to SLAG data, which include some of our range of interest, 
or lie close enough for extrapolation. We ignore the details of the low energy resonances 
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the relevant range, which are only approximately included in one of the parameterizations[|lj 
and not at all in the second |18[. We concentrate on the magnitude and shape of the back- 
ground from smoothly varying fits to the nucleon structure functions. Resonances should of 
course be added, but this will require more complete data separating the longitudinal and 
transverse contributions. 

The paper begins with a review of the basic formalism for (e, e') on complex targets, 
which defines the structure (or response) functions for nucleon or nuclear targets. In Section 
3 we formulate a simple impulse model for the response of a nuclear target based on a Fermi 
gas of bound nucleons. This follows earlier treatments of Fermi smearing, e.g., by Bodek 
and Ritchie, |20] but with some detailed differences in the prescription for going off shell. 
The treatment of binding effects continues in Section 4. The parametric fits to the nucleon 
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structure functions appear in Section 5. Finally, we present our estimates for the contribution 
of nucleon inelastic (internally excited) response functions in Section 6, and discuss these 
results in Section 7. 



2 Basic formalism 



In this section we summarize the basic formalism for (e, e') scattering from a nuclear target, 
which may be either a single nucleon or an A-hodj nucleus. The following sections give 
extensions for bound nucleons. In the Born (one-photon exchange) approximation one can 
write the differential cross section 





|k'| 


dn'dE' 


g4 |k| 



(2.1) 



where L^'^ is a lepton tensor describing incoming and outgoing plane-wave electron states, 
k and k' represent the initial and final electron four-momenta, and q = k — k' is the four- 
momentum transferred to the target via virtual photon exchange. 

Let p and p' represent the initial and final four- momenta of the target. In principle 
the target (response) tensor PF^j, is a function of all three variables p, p' and g, however, 
conservation of four- momentum can be used to eliminate any one of these; it is conventional 
to eliminate p' . The most general form of W^y satisfying Lorentz invariance, gauge invariance 



and parity can then be written [pi 



+ W2 



M 



M g2 



M 



p-qqv 

M g2 



(2.2) 



where the tensor behavior of W^^ under Lorentz transformations is described entirely by the 
quantities in square brackets. The scalar structure functions Wi and W2 depend only on 
scalar combinations of p and q (and the target mass M), from which one can form three 
scalar combinations: p^, q^ andp-g, or equivalently = — g^, u = p-q/M and W"^ = (p+q)^- 
For an on-shell target p"^ = so only two independent variables remain. A conventional 
choice for nucleons, which we adopt, is and u. 

For a nuclear target at rest in the laboratory frame, inserting (2^) into (|2.1|) leads to 



da. 



dQ'dE' dQ' 



(2.3) 



^l^,>,q)+ Q^ + tan2^)w^^^(.;,q) 

where the superscript A refers to the A-body nucleus. The longitudinal (W^) and transverse 
(W^) structure (or response) functions are defined 



Wiiu, q) = ^W.^iuj, q) - W.^iu, q) 



(2.4) 
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W^{u, q)=2x W^{uj, q) . (2.5) 

Equations (|2.3| )- (|2.5| ) have been written as functions of the lab variables u (energy transfer) 
and q (three-momentum transfer) with q'^ = {u, q), as is standard for analyzing the nuclear 
response. In these expressions the target mass M which appears in ( |2.2| ) has dropped out 
since Pq = M for an on-shell target at rest. Note also that in this case u = u. 

Since real photons are purely transverse, = at = 0, i.e., at u; = |q|. In order to 



analyze nuclear (e, e') data for the purpose of extracting two-body correlation functions |I5 



|T6[ the function must be evaluated over the entire range < u < |q|. Its behavior near 
= is therefore an important qualitative feature which must be maintained. 



3 Fermi smearing 

We are primarily interested in estimating the contribution of internal excited states of the 
nucleon to the inelastic nuclear response, compared to the nuclear response in which nu- 
cleons are not excited internally. For this comparison we need a model in which we can 
express the A-hodj target tensor W^i^{uj, q) in terms of W^^^^p, q), the corresponding 1-body 
tensor for constituent nucleons in the target of isospin projection cr. We also need quanti- 
tative information on the nucleon structure function, for both elastic and inelastic nucleon 
kinematics. 

Any realistic model which describes an A-hodj nucleus in terms of its constituent nucleons 
must include at least Fermi motion and nuclear binding effects. The simplest model is the 
plane- wave impulse approximation (PWIA) in which final state interactions are ignored. This 
should be adequate to provide a reasonable estimate of the inelastic nucleon background. 
Within the PWIA, West|2^ has shown that the A-hodj response tensor in the laboratory 
is given by [| 

where Ep = y/p'^ + M'^ and M is the nucleon mass. The momentum distribution n„{p) for 
nucleon species a is normalized to N„/2, where Np = Z and Nn = N, and the overall factor 
of 2 reflects a sum over both spin states. The factor Ep/M is required to preserve the phase 
space volume of d^p under Lorentz transformations PB| . 

For a free Fermi gas we would put (T^) on the right-hand side of ( ^.1| ), however, expression 



was derived for on-shell targets and is not strictly valid for bound nucleons. A complete 



^For elastic nucleons Pauli blocking in the final state should be included, and leads to an additional factor 



[1 - n^(p + q)] . See (Q and (g) 
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specification of for the general off-shell case is outside the scope of this paper because it 
requires knowledge of nucleon dynamics in the medium. We will simply assume that ( |2.2| ) is 
a valid starting point for an off-shell extension. Binding effects then enter PV^^ in two ways: 
1) through the initial-state nucleon energy po which appears explicitly in the square brackets 



in ( p.2|) , and 2) implicitly through the scalar functions W[ and We have already stated 
that on shell W[ and W2 are functions only two variables, which we take to be and 
u, and that for a bound nucleon W[ and W2 can in principle depend on three variables, 
e.g., Q^, V and W'^ . Since experimental data for W^^ is available only on shell, a minimal 
procedure is to assume that off shell and W2 are also functions of only two variables. 
In principle there exists considerable freedom to choose which two, however, in what follows 
we will show that and v is the natural choice in order to satisfy the requirement that 

= at g2 = 0. 

To procede we insert ( f^.2D into each side of (|3.1|). With the z-axis chosen along q and 
assuming spherical symmetry for n(p), equating tensor components leads to 



<(c.,q) = 2 E / 



(27r)3 Ep/M 



(3.2) 



1^2^(0;, q) = 2 I 



(Pp 



(27r)3Ep/M 



1 + 



M \q\u 



2 ,ly^2 



Inserting (U) and (jSj) into Q and (U) leads to 



m) q2 



W2^{Q',u). (3.3) 



d^p 



(27r)3 Ep/M 



M|q|z/ 



- H^2'^(g^//)-i^r(g^//) 



(3.4) 



d^p n^(p) 



(27r)3 Ep/M _ 

From (|3.4| ) it is easily verified that the condition = at = will be satisfied if 



(3.5) 



lim 



— w-2'(g^^/)-^^f(Q^z/) 



. 



(3.6) 



To understand the implications of (|3.6| ) we first observe that the factor v^jQ^ originates 

from the p ■ q terms in ( |2.2| ). Although the choice of variables on which W" and W2 depend 

off shell is in principle arbitrary, this work is based on the assumptions that 1) and 

W2 depend on only two variables, and 2) and W2 depend on the same two variables. 

^Technically speaking, the proton tensor Wjl^, is measured on shell by scattering from free protons, while 
the neutron tensor WJ^^, is determined by scattering from deuterons and extracting neutron contributions 
using theoretical arguments. 
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Consequently ( |3.6| ) will be satisfied only if and W2 are chosen to depend on Q'^ and v off 
shell. In contrast Bodek and Ritchie|2^, in a study of Fermi smearing at high Q^, assumed 
Wi and W2 to be functions of and W off shell, which is equivalent to choosing and v 
and evaluating v at the point vw = (W^^ — M"^ + Q'^)/2M. Since off shell i^vi/ 7^ z/ in general, 
this prescription will not satisfy (|3l6|) and implies 7^ at = 0. We therefore believe 
that their prescription is not useful when seeking an off-shell extrapolation valid at low . 
(Except for the region ~ 0, the effect of the two choices differs little.) 

When considering inelastic nucleons on shell it is convenentional to express Wi{Q'^,i>) 
in terms of H^2^(Q^,i^) and a third function R„[Q'^,v), which is proportional to the ratio 
Wl/Wt- 



1 + v^lQ" 



With the above assumptions for W[ and W2 , inserting 



:2D-(0) into 



and 



(3.7) 
I) gives 



(27r)3Ep/M 



Q-'V^ M\q\vJ \ujj l + R,{Q^,v) 



(27r)3 Ep/M 



From ( ^.71 ) we see that in order to satisfy ( ^.61 ) we must have R^^Q^, z/) = at 
Section 5 we provide a parameterized form for R^jiQ^ , v) which has this property. 



(3.9) 
0. In 



4 Spectator models of nuclear binding 

In order to evaluate ( p.8|) and ( p.9|) for bound nucleons we must specify a particular model 
which determines the off-shell kinematics. In conventional terminology, a "spectator" model 
specifies the energy p^ of a bound nucleon in terms of its three-momentum p, which is related 
to the excitation energy of the recoiling residual (A— l)-body target nucleus (spectator). This 
in turn determines the off-shell value of the variable v = p ■ q/M. As in the PWIA of (|3.1|) , 
a spectator model is based on the assumption that nuclear interactions enter only the initial 
state; the final excited nucleon state is free. In this section we discuss three such models for 
inelastic nucleon response. 

Bodek and Ritchie have investigated the effects of Fermi motion and nuclear binding 
beginning with ( |3.1|) . In that work they assumed a nucleon energy of the form 
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Po = Ma- v/p2 + Ml_i , (4.1) 

where Ma = A{M — E/A) is the mass of the A-hodj target nucleus and M^-i = Ma — M 
is the mass of the recoihng spectator nucleus after nucleon knockout. This form is a direct 
generalization of the expression for dissociation of the deuteron, as given by Atwood and 
West[^, and corresponds roughly to the "separation" energy, i.e., the energy required to 
remove a particle from the least bound state in the target. Compared to the free nucleon 
energy pq = Ep, which is bounded hj M < Ep < M + 34.5 MeV on the range < |p| < Pf, 
( [4.1| ) has relatively weak dependence on p, being bounded hj M < pq < M — 0.6 MeV 
over the same range.Q Expression ( |4.1| ) has two peculiar features when applied to a many- 
body system: 1) pq decreases with increasing p, and 2) pq includes essentially no binding 
effects, since po = M aX p = and for p 7^ the binding energy E/A enters only negligibly. 
Furthermore, because ( |4.1| ) corresponds to the separation energy it does not account for the 
possibility that the ejected nucleon was initially deeply bound. 

A simple modification of ([4.1|) is to account for the variation of po with p using a self- 
consistent potential model. We employ a relativistic mean field model based on quantum 
hadrodynamics (QHD) p4[, a relativistic quantum field theory of hadronic matter with 
Lorentz vector and scalar meson interactions. In this model the energy of a bound nucleon 
is given by 

Po = Vo + Vp2 + M*2 , (4.2) 

where M* = M + S* is the effective nucleon mass, 5* < represents an attractive scalar field 
and Vo > represents a repulsive, time-like vector field. The binding effects are much greater 
than in (^). The energy po in (|]|) is bounded by M - 76 MeV <Po< M -24 MeV over 
the range < |p| < pp, i.e., nucleons in the Fermi sea are bound by approximately 50 MeV. 
Furthermore, po has the expected behavior, i.e., increases with increasing p, and correctly 
accounts for the average binding energy for nucleons bound deeply in the Fermi sea. In this 
sense ( [4.2|) corresponds to the "removal" energy from occupied orbitals in the target. 

Simply using ( [4.2|) to evaluate z/ in ( |3.8|) and (|3.9| ) is not consistent, however, because 
( p. 21) was derived for a nucleon of mass M. In the mean field QHD theory, a bound nucleon 
acquires an effective mass M*, on which its Lorentz transformation properties are based. 
For a nucleon of mass M* ( |2.2| ) must be replaced by 

W^,{p,q) = W,^ 

•^The numbers quoted in this section are based on numerical parameters for ^^Fe given in Section 6. 



-9^,u + 



M* M* g2 



Pi 



P ■ q qu 

* „2 



M* M* q 



(4.3) 
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where p* = {E*, p) and = yp^ + M*^. (Notice that the vector field Vq does not enter the 
definition of p*, just as it does not enter a spinor Us{p) describing a Dirac plane- wave.) It is 
then natural to introduce the variable u* = p* ■ q/M* . Repeating the arguments of Section 
3, we assume that off-shell and W2 are given by the on-shell functions evaluated at 
and v* . This leads to the following expressions, which replace (|3.8|) and (|3.9|) 



d^p n^(p) [qV^ Pz 



1 + v*'IQ 



(4.4) 



W^^(^,q)=2x2 [ 



d^p n^(p) 
(27r)3E*/M* 



+ 



Px 



1 + v*VQ^ 



(4.5) 



These are the primary formulae used in our calculations of the inelastic nucleon background. 
These expressions are also used to compute the elastic nucleon response, i.e., quasielastic 
peak, for a bound Fermi gas. Numerical results, along with essential differences in interpre- 
tation, are given in Section 6. 



5 Nucleon structure functions 



In this section we describe the forms of the nucleon structure functions W^^ used in our 



calculations of the inelastic nucleon background. We rely on two parametric fits ||17| , [18 
which characterize (e, e') data obtained over a large kinematic range at SLAC, but which 
unfortunately do not cover all of the lower energy and momentum transfers available at 
CEBAF. We therefore must extrapolate these fits to this lower range, which introduces some 
uncertainty into the resulting estimates, as we shall see by comparing the results for different 
fits. 

The two experimental structure functions are often presented in terms of W2 and R^, as 
in ( ^.7| ), and given in terms of the invariant variables and z/, or and x = jlMv] 
the latter choice is of interest for describing scaling behavior at high x. Data for the proton 
functions comes directly from ^H(e, e')X, while the neutron functions must be extracted 
from ^H(e, e')X data using theoretical assumptions about the momentum distribution and 
binding effects in the deuteron. 



Bodek et. al. [l?! fit data for in the kinematic range 1.0 < g2 < 20 GeV^ and 



0.1 < a; < 0.77. They assume constant Rp = Rn = 0.18, which allows (|2.3|) to be expressed 



entirely in terms of VT^. Following a treatment developed by West|^, |2^, they use deuteron 



wavefunctions to account for Fermi smearing and binding effects in extracting the neutron 
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structure function. The parameterization for both proton and neutron structure functions 
is then given in the form 



W^{Q\ v) = B{Q\ W) W^{Q\ v) , (5.1) 

where B{Q^, W) is a modulating function for low energy excitations including threshold and 
resonant behavior, as a function of final invariant mass 

eliminate the resonant terms from this function since they do not appreciably modify the 
estimates, and since the value of R is expected to vary from resonance to resonance, contrary 
to the assumption of constant Rp = Rn- Without resonances, B{Q^^ W) vanishes below the 
vr-production threshold W ~ 1.07 GeV, increases with increasing W to a, second threshold 
near W ~ 1.74 GeV, and reaches unity for W ^ 2M; the functional form of B{Q'^,W) is 
given in Appendix A. The function W2 in ( p.lD is a smoothly varying function of and z/, 
and is also given in Appendix A. Its functional form was chosen to illustrate the behavior of 
VFf (Q^, z/) in a certain scaling variable (see Ref. 1171). 

Whitlowig fits in the kinematic range 0.6 < < 30 GeV^ and 0.06 <x< 0.9, 



using all SLAG data available at the time of publication. He first extracts R(j{Q^,x) using 
only data for which W > 2 GeV to exclude the resonance region. He finds Rp = Rd, and 
therefore -Rp = Rn, to within experimental errors (~ 5%), in good agreement with theoretical 
predictions based on QGD.|2^ However, noting that QGD predict ions [^, systematically 
underestimate the data, he provides a parameterization of the experimentally determined 
Rcr in the form 



RAQ^x) = ^ 



Ra{Q\x) + R,{Q\x) + R,{Q\x) 



for Q2 > 0.3 , (5.2) 



where Ra, Rb and Rc are given in Appendix B. 

For fixed x and ^ 1 GeV^, the experimental data indicates that R^j ~ 0.2 and 
decreases gradually as increases. For I GeV^, R^ decreases rapidly and tends 

to zero, as required to ensure Wl = at = 0. In contrast, if ( |5.2| ) is evaluated for 

< 0.3 GeV^ it will diverge as 0. To make reasonable estimates, therefore, we 

extrapolate Rcr{Q'^,x) down to = by assuming a simple analytic form which behaves 
properly at all Q^. A theoretically motivated formpHl is 



^'^(^''^) = r7T^' ioTQ'<0.3, (5.3) 

where Qr^x) and 6i?(x) are chosen such that ( |5.3| ) smoothly matches ( p.2[ ) at Q'^ = 0.3 GeV^. 
We have verified that the results in Section 6 are insensitive to the detailed form of this 
continuation, provided it matches ( |5.2| ) and tends smoothly to zero at = 0. 
The proton structure function W2 is parameterized in the form 
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X] 



1 + Ai(x) In 



A[x) 



+ Asfx) In^ 



A[x) 



(5.4) 



where explicit forms for the functions appearing here are given in Appendix C. The dominant 



behavior of 



is determined by the function F\ 



thr. 



X 



which is a power series in (1 — a;) 



and therefore tends to zero at the elastic nucleon threshold at x 
function is parameterized in terms of the ratio 



wi 



Piix) + P2ix) \nQ 



1. The neutron structure 



(5.5) 



'2 \''''2/s 

where the subscript S refers to Fermi smearing effects in the deuteron, which Whitlow points 
out are negligible for x < 0.5. In this region the data looks almost linear in x, although 
strong curvature is apparent for 0.5 < x < 0.7. This curvature is most likely an artifact 
of Fermi smearing in the deuteron. Whitlow gives two different parameterizations of the 
functions Pi{x) and P2{x), one essentially linear in x and one which shows some curvature 



corresponding to this Fermi motion effect. Bodek et. al. have extracted the non-smeared 
ratio W2/W2 and find linear x-dependence over their whole range of data. For the purpose 
of our estimates, therefore, we choose the more linear of Whitlow's two parameterizations 
even though it was fit only for x < 0.5. Explicit forms for the functions Pi{x) and P2{x) 
which are used in ( |5.5D are given in Appendix C. 

Before proceeding to our estimate of the inelastic nucleon background, it is useful to point 
out some general features of the fits to W2, since these largely determine the behavior of the 
nuclear response functions ( ^.4| ) and ( [4.5| ). Although data for W2 are shown in the original 
references, they are typically shown at higher energies and in terms of invariant variables, 
e.g., and x. Figures 1-3 show W2 and W2 in laboratory variables u and |q — , for free 
nucleons at rest and for three-momentum transfers |q| = 1, 2 and 4 GeV/c. The energy 
range shown is — M < uj < |q|, i.e., from the position of the elastic nucleon peak to the 
absolute upper limit attainable by electron scattering. 

We first focus on protons. At these relatively low momenta, the fit by Bodek et. al. (B) 
is significantly smaller than that by Whitlow (A12), with nearly a factor of two difference 
for |q| = 1 GeV/c. This difference decreases with increasing |q|, and by |q| = 4 GeV/c 
the two fits are nearly indistinguishable. At |q| = 2 GeV/c, B shows the second (higher) 
threshold behavior of (|A^), while at |q| = 1 GeV/c the values of the final invariant mass 
W are such that all u lie entirely below the second threshold. The behavior of B and A12 
at the endpoints differs qualitatively. At the lower limit, i.e., the elastic threshold, A12 
vanishes at u = i/q^ + M'^ — M (see discussion following ( |5.4| )), while B vanishes at the 
first inelastic threshold lo = ^/^~+W^ — M with W ^ 1.07 GeV/c. At the upper limit 
B vanishes precisely at a; = |q|, while the logarithmic behavior of A12 in must be cut 
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off to allow only positive values, causing W2 to vanish slightly prior to uj = |q|. Generally 
speaking, the fits to W2 follow similar trends, although with different numerical values. For 
both fits the neutron-to-proton ratio W2/W2 < 1. Compared to the B fit, this ratio in the 
A12 fit, i.e., that given by (|5.5|) is ~ 40% smaller at the lower limit and ~ 10% larger at 
the upper limit for |q| = 1 GeV/c, and is only ~ 25% smaller at the lower limit and ~ 5% 
larger at the upper limit for |q| = 2 GeV/c. In Fig. 1, the difference in shape between W2 
and W2 in the A12 fit is reflective of the fact that at lower |q| the slope of the ratio W2/W2 
is greater than that in the B fit. 



6 Numerical results 

In this section we present numerical results for the inelastic nucleon background obtained 
by using ( |4.4| ) and ( |4.5| ) with the fits to W2 and R^j given in Section 5. We assume a simple 
Fermi gas momentum distribution no-(p) = 9(pp — |p|) with pp = 0.257 fm^^, corresponding 
roughly to ^^Fe. The free nucleon mass M = 938.92 MeV. To account for nuclear binding 



we use a reduced mass M* /M = 0.648, corresponding to the mean field theory of Ref. p4 
for the same value of pp- This value of M* represents the smallest reduced mass which can 
reasonably be expected in a real nucleus, so that comparing this result with that for free 
nucleons demonstrates the range of results which can be expected in this model. In all cases 
where W and x must be replaced by the variables and u in the functions W2 and Ro^, as 
mentioned in Section 5, these replacements have been made in the on-shell functions before 
letting z/ — > z/*. 

In what follows the inelastic background will be shown compared to the corresponding 
quasielastic peak at the same three-momentum transfer, whose calculation for free nucleons is 
summarized in Appendix D. To include the effects of nuclear binding on the quasielastic peak 
we follow the theoretical treatment of Ref. |l^, which is described at the end of Appendix D. 
This leads to the standard expressions for quasielastic scattering from a Fermi gas composed 
of nucleons of mass M*, for which binding effects enter both initial and final nucleon states. 
Note that this differs from the spectator models described in Section 4, for which binding 
effects enter final states only indirectly, since there is assumed to be no final state interaction 
of the inelastically excited nucleon. Therefore, the comparison between the elastic and 
inelastic nucleon contributions computed here is given only to provide a reasonable estimate 
of their relative sizes and positions. 

Figure 4 shows the longitudinal (a) and transverse (b) nuclear response functions at 
three-momentum transfer |q| = 1 GeV/c. Thin curves are for free nucleons with mass M 
and thick curves are for bound nucleons with mass M*. As seen in Fig. 1, at this three- 
momentum the fit by Bodek et. al. (B) implies a background roughly half as large as that 
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implied by Whitlow's fit (A12). (Note that the energy range extends slightly lower than that 
in Fig. 1 to accommodate Fermi broadening.) As in Fig. 1, B vanishes at exactly uj = |q|, 
while A12 vanishes at slightly lower uo due to the logarithmic behavior of ( ^.4|) . Note that 
although the inelastic nucleon results derived from B and A12 differ, they are of the same 
magnitude as the quasielastic peak and displaced sufficiently that respective peaks appear 
at distinct energies. It therefore seems possible to effect a separation of elastic and inelastic 
nucleon contributions in this momentum range using these data fits, albeit with substantial 
errors reflecting the difference between the fits. 

The reduced mass M* affects the quasielastic peak much more than the inelastic back- 
ground. To understand this result we first note that Fermi motion is treated identically 
in the two calculations, and that in (|4.3| ) the effective mass always enters the coefficients 
of and W2 through the ratio Pq/M* = E^/M* ~ 1. Thus differences in sensitivity to 
the value of M* originate in the nucleon structure functions Wi and We have chosen 
to evaluate the inelastic nucleon structure functions in terms of and z/* ~ cj, hence the 
inelastic nucleon background is not strongly sensitive to the value of M*. It is in this sense 
that our background calculation includes binding effects only indirectly in the final state. In 
contrast, in order to match standard treatments the elastic nucleon structure functions ( pi[ ) 
and ( p2| ) have explicit factors of r* = Q'^/AM*'^, which scales quadratically with the reduced 
mass M*. Hence the size and location of the quasielastic peak is more sensitive to the value 
of M* than is the inelastic nucleon background. This is clearly an artifact of the particular 
model we have used for the quasielastic peak, and would change with any modification of the 
theory for either elastic or inelastic nucleons. Such changes are not of interest here, since our 
main goal is to provide a reasonable estimate the inelastic background and determine how 
reliably these parameterizations of the SLAG data can be extended to the CEBAF kinematic 
range. 

Figure 5 shows the longitudinal (a) and transverse (b) nuclear response functions at 
three-momentum transfer |q| = 2 GeV/c. The curve labeling is the same as in Fig. 4. The 
threshold behavior of B, which can be seen in Fig. 2, is not visible in Fig. 5 because of 
Fermi smearing. The most notable difference with Fig. 4 is that the inelastic background is 
much larger compared to the quasielastic peak. The scales of the plots show that this is due 
mainly to the rapid decrease of the dipole form factor (|D7|), and due to a lesser extent to 
a slight increase in the inelastic background. Thus separating elastic and inelastic nucleon 
contributions from experimental data will be much more challenging for |q| :> 2 GeV/c - 
an effect which has long been recognized but not quantitively investigated. However, this 
difficulty is partly compensated by the fact that, as noted in Fig. 2, the two inelastic fits 
B and A12 are in much closer agreement, and therefore the inelastic background can be 
evaluated with more confidence. Results for |q| = 4 GeV/c are not shown because by this 
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point the dipole form factor has totally suppressed the elastic peaks relative to the inelastic 
background. 



7 Summary and conclusions 

The study of nuclear structure by (e, e') reactions at a few GeV, the CEBAF kinematic range, 
requires the removal of the background resulting from inelastic excitation of single nucleons. 
In this paper we provide an estimate of this background based on two different parameteriza- 
tions of inelastic nucleon structure functions measured at SLAG. We then compare the sizes 
and positions of the inelastic nucleon background to those of the quasielastic peak, which 
represents the dominant contribution of elastic nucleons to inelastic nuclear excitation. We 
assess how confidently the available SLAG fits can be applied in the GEBAF kinematic range 
from the difference between the results for each fit over that kinematic range. 

We make a number of assumptions which allow us to relate inelastic nucleon excitations 
in the nuclear target to those of free nucleons; these all involve allowing only "minimal" 
effects of the nucleus on internal nucleon excitations. First we use the PWIA in ( |3.1|) , which 
introduces Fermi motion but does not include final state interations for the excited nucleon. 
Second we assume that the nucleon response tensor W^^ for a bound nucleon, given in ( |4.3| ) 
is of the same tensor form as that for a free nucleon, given in ( p.2|) , but modified for a nucleon 
with effective mass M* . Third we assume that the nucleon structure functions and W2 
for bound nucleons are equal to the free-nucleon structure functions evaluated at the same 
value of and dX u = u* (see discussion leading to ( [4.4|) and ([4.5|) ). Thus binding effects 



enter the enter final states in (|4.3| ) only indirectly, which is consistent with ignoring final 
state interactions in (|3.1| ). This rule for off-shell extrapolation has good behavior at = 0, 
as seen in (^]6|). 

Our calculation of the inelastic response functions uses parameters based on the structure 
of ^^Fe, but with a sharp Fermi distribution. There is little sensitivity of the inelastic nucleon 
background to the choice of M*, which enters only minimally. In contrast, our model for 
the quasielastic peak, which does have final state interactions, is more sensitive to M*, but 
this does not significantly affect our comparison of the relative size and position of the two 
contributions. We find that for |q| ^ 1 GeV/c the elastic and inelastic nucleon contributions 
are comparable in size, and can be separated, but in the inelastic nucleon response there is an 
uncertainty of roughly a factor of two resulting from disagreement between the SLAG fits in 
this kinematic range. For |q| ^ 2 GeV/c improved agreement between the two extrapolations 
raises the certainty of the background calculation, but the momentum dependence of the 
Sachs form factors quenches the elastic nucleon contribution with increasing |q|. 

Our estimates of the inelastic nucleon background show that with the presently avail- 
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able SLAC fits it is feasible to extract nuclear information from nuclear (e, e') data for 
|q| <^ 1 GeV/c. However, the present accuracy does not make this sufficiently useful, since 
many interesting effects (e.g., of correlations) make rather small contributions to the response 
functions. What is needed is more complete data and analysis on nucleon structure functions 
at lower Q^, with full separation of Wi and W2, i.e., R^. Clearly resonances are important 
in this kinematic region, and should be included in the inelastic nucleon background. In con- 
trast, for |q| ^ 2 GeV/c separating the inelastic nucleon background can be done with more 
confidence, but the momentum dependence of the Sachs form factors makes identification 
of the elastic nucleon response difficult. It seems that this situation can not be improved 
by making more precise measurements of the inelastic nucleon response functions, since the 
existing fits are already in good agreement in this kinematic range. For some purposes, the 
use of (e, e'p) and (e, e'n) experiments would reduce the problem of the background, but 
these reactions are usually less complete kinematically than (e, e'), and are more sensitive to 
final state interactions. These features limit their usefulness in sum rule studies. 
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Appendix A: Fit to W2{Q'^, i^) by Bodek et. al. 



In Ref. the proton and neutron structure functions are parameterized in the form 

g2 2Mz/ + ci 



2Mz/2 Q^ + C2 



7 

. n=3 



1 - 



2Mu + ci 



(Al) 



where Ci = +1.6421 and C2 = +0.3764. The remaining parameters c„o- are listed in Table 1. 





a = p 


a = n 




+0.2562 


+0.0640 


C4a 


+2.1785 


+0.2254 


C^u 


+0.8978 


+4.1062 


Cda 


-6.7162 


-7.0786 


Cla 


+3.7557 


+3.0549 



14 



Table 1. Numerical parameters for W2, taken from Bodek et. al. ||T7 
The modulating function B{Q^, W) is given by 



B{Q\ W) = B{W) [1 + (1 - B{W)) [h + h{x - hf)_ 



where x = Q'^/{W^ - + Q^) and B{W) is defined 



B{W) = 9{W-bi) 62 



+ 9{W - 64) (1 - &2) 



1 _ Q-h{w^ 



The coefficients 6j are listed in Table 2. 



= +1.0741 


62 = -0.7553 


63 = +3.3506 


64 = +1.7447 


65 = +3.5102 


6g = -0.5999 


67 = +4.7616 


bs = +0.4117 



Table 2. Numerical parameters for B{Q'^, W), supplied by Bodek [p5|] 

Appendix B: Fit to R{Q'^,x) by Whitlow 



The functions appearing in ( ^.2|) are given by 

b2 



Rc{Q\x) 



ln(gV0.04) 
bi 



1/4 



ln(gV0.04) 

Cl 



e(g^x) + -4 + 



(g4 + 0.3^ 



in(gvo.o4) 



Q{Q\x) + 



C2 



(g2-5(i-x)5) +C 



where 



e(g^x) = 1 + 12 



0.1252 



.g2 + iyva;2 + 0.1252 

and the coefficients appearing in (pTD-(p3|) are listed in Table 3. 



ai = +0.0672 


a2 = +0.4671 


as = +1.8979 


61 = +0.0635 


62 = +0.5747 


63 = -0.3534 


Cl = +0.0599 


C2 = +0.5088 


C3 = +2.1081 
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Table 3. Numerical parameters for R{Q^, x), taken from Whitlow p8[ . 



Appendix C: Fit to W^iQ^^x) by Whitlow 



The functions appearing in ( ^.41) are given by 



Ft{Q',x) = Y.d, (1-x) 

i=l 



i+2 



A2(x) 



1=0 

Etodi+ex' ifQ^<Aix) 
otherwise , 



A{x) = 1.22 e' 



3.2x 



where the coefficients di are listed in Table 4. 



di = +1.417 


d2 = -0.108 


rfg = +1.486 


= -5.979 


dr, = +3.524 


de = -0.011 


dj = -0.619 


ds = +1.385 


dg = +0.270 


dio = -2.179 


du = +4.722 


d^2 = -4.363 



Table 4. Parameters for Whitlow's A12 fit to W2{Q'^,x), taken from Ref. |T8 
The functions Pi{x) and ^2(2;) appearing in ( ^.5|) are given by 



Pi (a;) ~ 0.9498 - 0.9706 x + 0.3102 x^ 



P2{x) ^ -0.0146 

where the numerical values are taken from Ref. |l^ . 



Appendix D: Elastic nucleon structure functions 

The elastic structure functions for free nucleons may be written 
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where (^^^(Q^) and ^^/^(Q^) are the Sachs electric and magnetic form factors, and r = 
QV4M2. Inserting (0) and O into Q and Q leads to 

W^(uj, q) = 2 X 2 / ^ /„(p) B^^Si^S(^ + Bp - Bp+J , (D4) 

where 



l + r 



-{Ep + Ep+q 



,2 , ^llaiQ' 



l + r 



rjq 



(D5) 



^n(p,q;Q' 



l + r 



-(2 p.)' + 



l + r 



r(2 pj2^(l + r)g^ 



(D6) 



The factor /o-(p) = '^(t(p) 1 — '^(t(p + q) in (P3|) and (|D4|) provides Pauli blocking for final 
states, and must be included when starting from (|3.1| ) (see footnote 1). Expressions ( P3|) 



and (D4) can be evaluated numerically by rewriting the (5- function as a 6'-function which 
restricts the polar angle = cos^"*^ [p ■ q/|p||q|]. 

In numerical calculations we take the Sachs form factors to be 



GeAQ^) = (l + gV0.71GeV2)-2, (D7) 
GmM) = (1 + t^p)GEM) > (D8) 

GMn{Q^) = K^GEpiQ^) , (D9) 
GEniQ^) = 0. (DIO) 

where Kp = +1.79 and = — 1.91 are the proton and neutron anomalous magnetic moments, 
respectively. 

For nucleons of mass M* the above expressions must be modified. This is accomplished in 
part by letting Ep Ep, -Ep+q — > E*^^ and r — > r* in expressions (P3|)-(P6|). While these 
changes must be made to account for the modified kinematics of nucleons with a reduced 
mass M*, there is some ambiguity in the treatment of the Sachs form factors which are fit 



only to free nucleon data. In this work we adopt Model G of Ref. [IC], in which we assume 



that the Sachs electric and magnetic form factors are unmodified in the nuclear medium. 
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Figure Captions 



FIG. 1. Fits to nucleon structure functions W2 and W2 in laboratory variables for three- 
momentum transfer |q| = 1 GeV/c. Fits by Whitlow (A12) for the proton (solid) and neutron 
(dashed), and fits by Bodek et. al. (B) for the proton (dot-dashed) and neutron (dotted). 

FIG. 2. Same as Fig. 1 except for |q| = 2 GeV/c. 

FIG. 3. Same as Fig. 1 except for |q| = 4 GeV/c. 

FIG. 4. Nuclear response functions (per nucleon) for three-momentum transfer |q| = 
1 GeV/c: (a) longitudinal W^{u;,q)/A, and (b) transverse W'^{u;,q)/A. Inelastic nucleon 
background based on fits by Whiltow (A12) (solid curves) and by Bodek et. al. (dashed 
curves) are shown, along with the quasielastic peak (dotted curves). Thin lines are for free 
nucleons with M*/M — 1, and thick lines are for interacting nucleons with M*/M — 0.648. 

FIG. 5. Same as Fig. 4 except for |q| = 2 GeV/c. The inserts show enlarged views of the 
quasielastic peak region. 
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